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■ Abstract 

o ■ 

04 ' Positive frequency Wightman function and vacuum expectation value of the energy- 

momentum tensor are computed for a massive scalar field with general curvature coupling 
parameter subject to Robin boundary conditions on two parallel plates located on (D + 1)- 
i dimensional AdS background. The general case of different Robin coefficients on separate 

^\ ' plates is considered. The mode summation method is used with a combination of a variant 

. of the generalized Abel-Plana formula for the series over zeros of combinations of cylinder 

' functions. This allows us to extract manifestly the parts due to the AdS spacetime without 

. boundaries and boundary induced parts. The asymptotic behavior of the vacuum densities 

«) ' near the plates and at large distances is investigated. The vacuum forces acting on the 

, boundaries are presented as a sum of the self-action and interaction forces. The first one 

"•pi ' contains well-known surface divergences and needs further regularization. The interaction 

ji. . forces between the plates are attractive for Dirichlet scalar. We show that there is a region 

^ ' in the space of parameters defining the boundary conditions in which the interaction forces 

, are repulsive for small distances and attractive for large distances. An application to the 

^ ' Randall-Sundrum braneworld with arbitrary mass terms on the branes is discussed. 

, 

H ! 1 Introduction 

, 

Anti-de Sitter (AdS) spacetime is one of the simplest and most interesting spacetimes allowed by 
general relativity. Quantum field theory in this background has been discussed by several authors 
(see, for instance, Refs. [1]-[17]). Much of early interest to AdS spacetime was motivated by the 
questions of principle related to the quantization of fields propagating on curved backgrounds. 
The importance of this theoretical work increased when it was realized that AdS spacetime 
emerges as a stable ground state solution in extended supergravity and Kaluza-Klein models 
and in string theories. The appearance of the AdS/CFT correspondence and braneworld mod- 
els of Randall-Sundrum type has revived interest in this subject considerably. The AdS/CFT 
correspondence (for a review see [18]) represents a realization of the holographic principle and 
relates string theories or supergravity in the bulk of AdS with a conformal field theory living on 
its boundary. It has many interesting formal and physical facets and provides a powerful tool 
to investigate gauge field theories, in particular QCD. Recently it has been suggested that the 
introduction of compactified extra spatial dimensions may provide a solution to the hierarchy 
problem between the gravitational and electroweak mass scales [19, 20, 21]. The main idea to 
resolve the large hierarchy is that the small coupling of four dimensional gravity is generated 
by the large physical volume of extra dimensions. These theories provide a novel setting for 
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discussing phenomenological and cosmological issues related to extra dimensions. The model 
introduced by Randall and Sundrum is particularly attractive. Their background solution con- 
sists of two parallel flat branes, one with positive tension and another with negative tension 
embedded in a five dimensional AdS bulk [20]. The fifth coordinate is compactified on 5*^/^2, 
and the branes are on the two fixed points. It is assumed that all matter fields are confined on 
the branes and only the gravity propagates freely in the five dimensional bulk. In this model, 
the hierarchy problem is solved if the distance between the branes is about 40 times the AdS 
radius and we live on the negative tension brane. More recently, scenarios with additional bulk 
fields have been considered [22, 23, 24, 25, 26]. 

In the scenario presented in [20] the distance between the branes is associated with the 
vacuum expectation value of a massless scalar field, called the radion. This modulus field has 
zero potential and consequently the distance is not determined by the dynamics of the model. 
For this scenario to be relevant, it is necessary to find a mechanism for generating a potential 
to stabilize the distance between the branes. Classical stabilization forces due to the non-trivial 
background configurations of a scalar field along an extra dimension were first discussed by Gell- 
Mann and Zwicbach [27]. With the revived interest in extra dimensions and braneworlds, as 
modified version of this mechanism, which exploits a classical force due to a bulk scalar field with 
different interactions with the branes, received significant attention [28] (and references therein). 
Another possibility for the stabilization mechanism arises due to the vacuum force generated 
by the quantum fiuctuations about a constant background of a bulk field. The braneworld 
corresponds to a manifold with dynamical boundaries and all fields which propagate in the 
bulk will give Casimir-type contributions to the vacuum energy (for the Casimir effect see Refs. 
[29]-[34]), and as a result to the vacuum forces acting on the branes. In dependence of the 
type of a field and boundary conditions imposed, these forces can either stabilize or destabilize 
the braneworld. In addition, the Casimir energy gives a contribution to both the brane and 
bulk cosmological constants and, hence, has to be taken into account in the self-consistent 
formulation of the braneworld dynamics. Motivated by these, the role of quantum effects in 
braneworld scenarios has received some recent attention. For a conformally coupled scalar this 
effect was initially studied in Ref. [35] in the context of M-theory, and subsequently in Refs. 
[36]-[49] for a background Randall-Sundrum geometry (for the related heat kernel expansions 
see Refs. [50, 51, 52]). The models with dS branes are considered as well [48, 53, 54, 55]. 
For a conformally coupled bulk scalar the cosmological backreaction of the Casimir energy is 
investigated in Refs. [35, 48, 53, 56, 57, 58]. 

Investigation of local physical characteristics in the Casimir effect, such as expectation value 
of the energy-momentum tensor, is of considerable interest. In addition to describing the physical 
structure of the quantum field at a given point, the energy-momentum tensor acts as the source 
in the Einstein equations and therefore plays an important role in modelling a self-consistent dy- 
namics involving the gravitational field. In this paper we will study the vacuum expectation value 
of the energy-momentum tensor of a scalar field with arbitrary curvature coupling parameter 
obeying Robin boundary conditions on two parallel plates in (Z)-|-l)-dimensional AdS spacetime. 
The general case is considered when the constants in the Robin boundary conditions are different 
for separate plates. Robin type conditions are an extension of Dirichlet and Neumann boundary 
conditions and appear in a variety of situations, including the considerations of vacuum effects 
for a confined charged scalar field in external fields [59], spinor and gauge field theories, quantum 
gravity and supergravity [60, 61]. Robin conditions can be made conformally invariant, while 
purely-Neumann conditions cannot. Thus, Robin-type conditions are needed when one deals 
with conformally invariant theories in the presence of boundaries and wishes to preserve this 
invariance. It is interesting to note that the quantum scalar field satisfying the Robin condition 
on the boundary of the cavity violates the Bekenstein's entropy-to-energy bound near certain 
points in the space of the parameter defining the boundary condition [62] . The Robin boundary 
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conditions are an extension of those imposed on perfectly conducting boundaries and may, in 
some geometries, be useful for depicting the finite penetration of the field into the boundary with 
the 'skin-depth' parameter related to the Robin coefficient. Mixed boundary conditions natu- 
rally arise for scalar and fermion bulk fields in the Randall-Sundrum model [22, 43]. To obtain 
the expectation values of the energy- momentum tensor, we first construct the positive frequency 
Wightman function. The application of the generalized Abel-Plana formula to the correspond- 
ing mode sum allows us to extract manifestly the boundary-free AdS part. The expressions for 
the boundary induced vacuum expectation values of the energy-momentum tensor are obtained 
by applying on the subtracted part a certain second order differential operator and taking the 
coincidence limit. Note that the Wightman function is also important in consideration of the 
response of particle detectors at a given state of motion (see, for instance, [63]). 

We have organized the paper as follows. In the next section we consider the vacuum in the 
region between two parallel plates. The corresponding positive frequency Wightman function 
is evaluated by using the generalized Abel-Plana summation formula for the series over zeros 
of a combination of cylinder functions. This allows us to present the boundary induced part in 
terms of integrals with exponential convergence for the points away the boundaries. In Section 
3 we consider the vacuum expectation value of the energy-momentum tensor for the case of a 
single plate. Both regions on the right and on the left from the plate are investigated. Various 
limiting cases are discussed. Section 4 is devoted to the vacuum energy- momentum tensor for 
the geometry of two parallel plates. The corresponding vacuum expectation values are presented 
in the form of the sum of single plates and 'interference' parts. The latter is finite everywhere 
including the points on the boundaries. The interaction forces between the plates are discussed 
as well. In Section 5 we show that the vacuum expectation value for the energy-momentum 
tensor of a bulk scalar in the Randall-Sundrum braneworld scenario is obtained from our results 
as a special case. The last section contains a summary of the work. 



2 Wightman function 

Consider a scalar field (p{x) on background of a {D + l)-dimensional plane-symmetric spacetime 
with the line element 

ds^ = gikdxUx'' = e-'^'^^y^'n^^dx^'dx'' - dy"^, (2.1) 

and with rj^y = diag(l, —1, . . . , —1) being the metric for the D-dimensional Minkowski space- 
time. Here and below i,k = 0, 1, . . . , D, and /x, i/ = 0, 1, . . . , D — 1. The corresponding field 
equation has the form 

(^'''ViVfe + rr? + CR) ip{x) = 0, (2.2) 

where the symbol Vj is the operator for the covariant derivative associated with the metric gik, 
R is the corresponding Ricci scalar, and C, is the curvature coupling parameter. For minimally 
and conformally coupled scalars one has C = and C, = C,c = {D — 1)/4Z) correspondingly. Note 
that by making a coordinate transformation 

z = j e^^^Uy, (2.3) 

metric (2.1) is written in a conformally-flat form ds^ = e^'^'^ rjikdx"^ dx^ . 

Below we will study quantum vacuum effects brought about by the presence of parallel 
infinite plane boundaries, located at y = a and y = b, a < b, with mixed boundary conditions 

(4 + Bydy) ^{x) = 0, y = a,b, (2.4) 
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and constant coefficients Ay, By. The presence of boundaries modifies the spectrum for the 
zero-point fluctuations of the scalar field under consideration. This leads to the modification of 

the vacuum expectation values (VEVs) of physical quantities to compared with the case without 
boundaries. In particular, vacuum forces arise acting on the boundaries. This is well known 
Casimir effect. As a first stage, in this section we will consider the positive frequency Wightman 
function defined as the expectation value 

G+{x,x') = {0\ip{xMx'm, (2.5) 

where |0) is the amplitude for the vacuum state. In the next section wc use this function to 
evaluate the VEVs for the energy-momentum tensor. Note that the Wightman function also 
determines the response of particle detectors in a given state of motion. By expanding the field 
operator over eigenfunctions and using the commutation relations one can see that 

G+{x,x') = Y,^a{x)<pU^'), (2.6) 

a 

where a denotes a set of quantum numbers, and {ipaix)} is a complete set of solutions to the 
field equation (2.2) satisfying boundary conditions (2.4). 

On the base of the plane symmetry of the problem under consideration the corresponding 
eigenfunctions can be presented in the form 

faix') = M^nUy), (2.7) 

where (pkix^) are the standard Minkowskian modes on the hyperplane parallel to the plates: 

0k(x'^) = ^====, k^ = {u;,k), a;=V^^T<, k = \k\. (2.8) 

Here the separation constants are determined by boundary conditions (2.4) and will be given 
below. Substituting eigenfunctions (2.7) into the field equation (2.2) for the function fn{y) one 
obtains the following equation 



dy 



(e-^-^) + K + CR) fn = my/n. (2.9) 



For the AdS geometry one has a{y) = k^y, z = e'^^y^ /ku, and R = —D{D + 1)A;|,, where the 
AdS curvature radius is given hy l/ko- In this case the solution to equation (2.9) for the region 
a < y < 5 is 

Uy) = cne^'^/' [Mmnz) + KY,{mnz)] , (2.10) 
where Jv{x), Yy{x) are the Bessel and Neumann functions, and 

(L>/2)2 - D{D + 1)C + m'^/kl. (2.11) 

Here we will assume values of the curvature coupling parameter for which v is real. For imaginary 
u the ground state becomes unstable [4] . Note that for a conformally coupled massless scalar one 
has = 1/2 and the cylinder functions in Eq. (2.10) are expressed via the elementary functions. 
From the boundary condition on y = a one finds 

6. = -4^^^, z,=e''^^^/kn, j = a,b, (2.12) 

'{mnZa) 
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where we use the notation 

F^^\x) = AjF{x) + BjxF'ix), Aj = Aj + BjkDD/2, Bj = Bjkn, j = a,b (2.13) 

for a given function F{x). Note that for Neumann scalar one has Aj = BjD/2 and 

F^^\x) = = x'-^'V/^Fix)]', Aj = 0. (2.14) 

Prom the boundary condition on the plane y = b we receive that the eigenvalues m„ have to be 
solutions to the equation 

Cf{Zb/Za, rUnZa) = (m„Za)yj^) (m^Zfc) - yj«) {rUnZaVi^^ (m^Zfe) = 0. (2.15) 

We denote hy z = "yv^m n = 1, 2, . . ., the zeros of the function C^*(r7, z) in the right half-plane of 
the complex variable z, arranged in the ascending order, 7,^^„ < ^v,n+i- The eigenvalues for m„ 
are related to these zeros as 

rrin = kDl^^ne'"^"-'^ = -iv,nlZa- (2.16) 
The coefficient c„ in Eq. (2.10) is determined from the orthonormality condition 

dye^^-''^ jn{y)fn'{3)) = Snn' (2-17) 



/ 

•J a 



and is equal to 



2 

(^l=7^Yu^"^\u)TfiV,u), U = 7.,n, V = Zb/Za, (2.18) 



where we have introduced the notation 

Note that, as we consider the quantization in the region between the branes, Za < z < Z},, the 
modes defined by (2.10) are normalizable for all real values of u from Eq. (2.11). 

Substituting the eigenfunctions (2.7) into the mode sum (2.6), for the expectation value of 
the field product one finds 



(0|^(x)v.(x')|0) = ^_^3 , / dke^'^(---')^/i.(7.,„)T;^,?,7.,n), (2.20) 



k'^-\zz')D'^ 

n=l 

where x = {x^,x^,... ,x^~^) represents the coordinates in {D — l)-hyperplane parallel to the 
plates and 

Kiu) = -=====— g^{u,uz/za)g^{u,uz'/za). (2.21) 

V" /^S + k^ 

Here and below we use the notation 

g,{u,v) = Mv)Yj''Hu) - jW(«)y,(i;). (2.22) 

To sum over n we will use the summation formula derived in Refs. [64, 65] by making use of 
the generalized Abel-Plana formula [66, 64]. For a function h{u) analytic in the right half-plane 
Re n > this formula has the form 



h{u)H^^^\'nu) 

[C-\ii,u)0''\u) ^ 
TT K'i'\r]x) [/i(xe^^/2) + h{xe-^'/^)\ dx 



yEM7..)T.(.,7..) = I j(»)^(,) + y(^)^(,)-2^^^"- 



4 Jo Ri"-^ (x) Ri"-^ {x)li^^ ivx) - Ri^^ (?7x)li"^ (a;) ^ 



(2.23) 
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where Ii,{u) and K^{u) are the Bessel modified functions. Formula (2.23) is vaUd for functions 
h{u) satisfying the conditions 

\h{u)\ < ei{x)e^^\y\ \u\ ^ oo, u = x + iy, (2.24) 

and 

h{ue^') = -h{u) + o{u-^), u ^ 0, (2.25) 

where ci < 2(?7 — 1), x^''^<'°~'^£i(x) — > for x +oo. Using the asymptotic formulae for the 
Bessel functions for large arguments when v is fixed (see, e.g., [67]), we can see that for the 
function hi,{u) from Eq. (2.21) the condition (2.24) is satisfied \i z + z' + \t — t'\ < 2zi,. In 
particular, this is the case in the coincidence limit t = t' for the region under consideration, 
Za < z,z' < zi,. As for \u\ < k one has hy{ue'^^) = —h^{u), the condition (2.25) is also satisfied 
for the function hu{u). Note that hu{u) ~ n^^^*:" for n ^ and the residue term on the right of 
formula (2.23) vanishes. Applying to the sum over n in Eq. (2.20) formula (2.23), one obtains 



{0\ip{xMx')\0) = ^ ^ ^ 



2D^D-l 



^j^gjk(x-x') 



hu{u)du 



Ji''^\u) + Yr\u) 



(a)2, 



(2.26) 



TT 



duu 



i^auiuZa,UZb) 



A;2 



G^") {uZa, Uz)G(f) {uZa, Uz') COSh [V^^ - P(t - t') 



where we have introduced notations 

G^Hu, v) = Uiv)Kjj\u) - fJ\u)K,{v), j = a, b, 



K^\v)/K^^\u) 



Ki''\u)Il-'[v] 



(2.27) 
(2.28) 



(the function with j = h will be used below) . Note that we have assumed values of the coefficients 
Aa and for which all zeros for Eq. (2.15) are real and have omitted the residue terms in the 
original formula in Refs. [64, 65]. In the following we will consider this case only. 
To simplify the first term in the figure braces in Eq. (2.26), let us use the relation 

gy{u, uz/za)gu{u, uz' jza) 



{a)2. 



J^{uz/Za)Jv{uz' /Za) 



(2.29) 



2 ■' — ' u^t^^'UnA 

with Hjf\z), (5 = 1,2, being the Hankcl functions. Substituting this into the first integral in 
the figure braces of Eq. (2.26) we rotate the integration contour over u by the angle 7r/2 for 
/3 = 1 and by the angle — 7r/2 for j3 = 2. Under the condition z + z' — \t — t'\ > 2za, the integrals 
over the arcs of the circle with large radius vanish. The integrals over {0,ikza) and (0, —ikza) 
cancel out and after introducing the Bessel modified functions one obtains 



hiy{u.)du 



4»)'(n) + yr^(t.) 



(a)2. 







2£| r 

TT Jk 



duu , — Juiuz) Ji, (uz' ) 

li''\uza) K,{uz)K, [UZ 



(2.30) 



du u _ , . . 

Ik Kr\uZa) VV^^ 



cosh 



Vu^-k^{t-t'] 



Substituting this into formula (2.26), the Wightman function can be presented in the form 



{0\cp{x)ip{x')\0) = {Os\cp{x)ip{x')\Os) + ((^(x)^(x'))(") 

a I UZii 



k%-\zz')^l^ 



^]^gik(x-x') ^2.31) 



Ik 



duu- 



G^f ) {uZa, uz)G^f^ {uZa, uz') cosh \f^^fl^{t - t') 
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Here the term 



{Os\v{xMx')\Os) = Vvr^-i y ^''^"^^^'^ Vo V^?^ J.(uz)J,K) (2.32) 

does not depend on the boundary conditions and is the Wightman function for the AdS space 
without boundaries. The second term on the right of Eq. (2.31), 

h.D-l( /\D/2 f 

Mx)v,(x'))(") = - V Cz. ydke--(— ')x 

\/v? -k^it-t') , (2.33) 



ll '{uza) Kt,{uz)K^{uz') 
duu — , — - cosh 



does not depend on the parameters of the boundary at ^ = and is induced in the region 
z > Za hy a single boundary at z = Za when the boundary z = is absent. 

Note that expression (2.32) for the boundary-free Wightman function can also be written in 
the form 

POO 

{Os\ipix)ip{x')\Os) = k^-\zz'f^^ / dmmGi^ (xf',x'f';m)Mmz)Mmz'), (2.34) 

Jo 

where G\j^{x^ ,x'^;m) is the Wightman function for a scalar field with mass m in the D- 
dimensional Minkowski spacetime M^. The right hand side in Eq. (2.34) can be further sim- 
plified by using the expression 



21 



^^^^ ' '^"(27r)^/2 [(x-x02-(t-i'-i£)2]P-2)/4 ' ^^"^^^ 

with e > 0. Substituting this into formula (2.34) and making use of the integration formula 
from [68], one obtains 

(05|^(x)^(.')|05) = (2vr)P+'i)/2 ^-^^-^^^^^^(^^ - ^)^'-''^"Qt-.]l'\-) (2-36) 

k^-^T{v + D/2)v-''-^l^ ^ (D + 2V + 2 D + 2u ^ \\ 

~ 2-Pi -. , — — + , 



where Q^{v) is the associated Legendre function of the second kind, 2F\{a^h\c\u) is the hyper- 
geometric function (see, for instance, [67]), and 

^^^^(.-/f + (x-xy-(,-,--.sf^ (2.37) 

2zz' 

It can be checked that in the limit fc^j — > from (2.32) the expression for the (Z) + l)-dimensional 
Minkowski Wightman function is obtained. To see this note that in this limit v ~ m/kD oo, 
k^z ~ 1 -|- koy and, hence, as it follows from (2.37), one has v ^ 1 + k'j^w'^ /2, where w"^ = 
{y ~ y')"^ + (^ ~ ^')^ — (t — t' — ie)2. Using the integral representations for the functions Q'^{v) 
and Kn{v) (formulae 8.8.2 and 9.6.23 in Ref. [67]), it can be seen that 



lim 



kl^e-'^-'iv' - l)-"/'Q^_i/2(^)l = ^K^imw). (2.38) 



For a = (D — l)/2 the expression on the right of this formula coincides with the {D + 1)- 
dimensional Minkowskian Wightman function up to the coefficient (27r)~^^~^^^^'^ . Hence, in 
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combination with the first formula in (2.36), relation (2.38) proves our statement. Alternatively, 
we can use formula (2.32), replacing the Bessel functions by their asymptotic expressions for 
large values of the order. 

In the coincidence limit x = x' expression (2.36) is finite for D <1 and for the VEV of the 
field square one obtains 

(o.i.-wio.>^ ^J;;,/'ff>;;'°^^^+;' . (2.39) 

This quantity is independent on the spacetime point, which is a direct consequence of the 
maximal symmetry of the AdS bulk. Formula (2.39) coincides with the result of Refs. [10, 14] 
(a typo of [10] is corrected in [14]) obtained from the Feynman propagator in the coincidence 
limit (for the zcta function based calculations see [12, 16]). The expression on the right of 
Eq. (2.39) is analytic in the complex D-plane apart from simple poles coming from the gamma 
function in the nominator. Hence, it can be extended throughout the whole complex plane. For 
even D this expression is finite and according to the dimensional regularization procedure [63] 
can be taken as the regularized value for the field square. Formula (2.39) can be also obtained 
directly from (2.32) firstly integrating over k and by making use of the integration formula 

I iD) = r dxx^-^jHx) - r(i^)r(z^/2 + .) 

By using the identity 



K^{uz)K^{uz') + Q.auiuZa,UZh)Gl"-^ (uZa,Uz)Gl^\uZa,Uz') 



^^^r-^^^^Ii,{uz)Iuiuz') + Qb„(uZa,UZb)G^j^\uZb,Uz)G^^\uZb,Uz'), 

iriuzb) 



(2.41) 



with 



nbu{u,v) = -^^ £^(^)/i%) (2.42) 

it can be seen that the Wightman function in the region Za < z < Zb can also be presented in 
the equivalent form 

{0\^{xMx')\0) = {OsMxMx')\Os) + ifixMx'))^'^ - ""^d^iJd J ^ke''-^-- ) (2.43) 
X r duu^^^^^^t^G^J'\uzb, uz)G^^^ {uzb, uz') cosh \^A^?^{t - t' 

where 



oo 



Ki^\uzb) Iv{uz)Ii,{uz') 



duu — , — cosh 



\lu^-k'^(t - t')l (2.44) 



is the boundary part induced in the region z < by a single plate at z = when the plate 
z = Za is absent. Note that in the formulae given above the integration over angular part can 
be done by using the formula 

/ Ae*(--V(*) = (2.)(^-)/^ 6k k-Mk) ■ (2.45) 
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for a given function F{k). Combining two forms, formulae (2.31) and (2.43), we see that the 
expressions for the Wightman function in the region Za < z < Zf, is symmetric under the 
interchange a b and I^, ^ K^. Note that the expression for the Wightman function is not 
symmetric with respect to the interchange of the plate indices. The reason for this is that, 
though the background AdS spacetime is homogeneous, the boundaries have nonzero extrinsic 
curvature tensors and two sides of the boundaries are not equivalent. In particular, for the 
geometry of a single brane the VEVs differ for the regions on the left and on the right of the 
brane. Here the situation is similar to that for the case of a spherical shell on background of the 
Minkowski spacetime. 



3 Casimir densities for a single plate 

In this section we will consider the VEV of the energy-momentum tensor for a scalar field in 
the case of a single plate located at z = Za- As it has been shown in the previous section the 
Wightman function for this geometry is presented in the form 

{0\^{xMx')\0) = {Os\v{xMx')\Os) + {^{xMx'))^-\ (3.1) 

where |0) is the amplitude for the corresponding vacuum state. The boundary induced part 
{(p{x)(p{x'))^"'^ is given by formula (2.33) in the region z > Za and by formula (2.44) with 
replacement z^ — ^ Za in the region z < Za- For points away from the plate this part is finite 
in the coincidence limit and in the corresponding formulae for the Wightman function wc can 
directly put x = x' . Introducing a new integration variable v = \/u^ — fe^, transforming to the 
polar coordinates in the plane {v, k) and integrating over angular part, the following formula 
can be derived 



^ .D-2r,^. uf{u) _0Fr(i^) r.....D-i 



dkk"^-^ / du = / duu^'-'fiu). (3.2) 

Jk V^?^ 2T{D/2) Jo ■'^ ' ^ ^ 

By using this formula and Eq. (2.33), the boundary induced VEV for the field square in the 
region z > Za is presented in the form 

The corresponding formula in the region z < Za is obtained from Eq. (2.44) by a similar way 
and differs from Eq. (3.3) by replacements Ii, ^ K^,. 

By using the field equation, the expression for the metric energy-momentum tensor of a 
scalar field can be presented in the form 



(^C - \^ gik^N^ - (ViVk - (Rik 



^\ (3.4) 



By virtue of this, for the VEV of the energy-momentum tensor we have 



(o|rifc(x)|o) = hm did'k{Q\ip{xMx'm + 



{o\^\xm. 



(3.5) 

This corresponds to the point-splitting regularization technique. VEV (3.5) can be evaluated 

by substituting expression (3.1) for the positive frequency Wightman function and VEV of the 
field square into Eq. (3.5). First of all we will consider the region z > Za- The vacuum 
energy-momentum tensor is diagonal and can be presented in the form 

{0\Tt\0) = {Os\Tl^\Os) + {Th^^K (3.6) 
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where 

is the VEV for the energy-momentum tensor in the AdS background without boundaries, and 
the term 

is induced by a single boundary at z = Za- For a given function g(y) the functions [g{v)] in 
formula (3.8) are defined as 

F^'Haiv)] = {\- 2C) [v'g'Hv) +(^D + vg{v)g'{v)+ (3.9) 

+{v'' + u'')g\v)]+{z^e-v'')9\v) 
F»[5(^)] = F(''^[g{v)] + (^v'-^^y\v), i = l,...,D-l, (3.10) 

+^ + + 2CI)(D + 1) - dV2] /(^), (3.11) 

and the expressions for the functions /^^'^[giv]] are obtained from those for F^'^^[g{v)] by the 
replacement v iv. Note that the boundary-induced part (3.8) is finite for the points away the 
brane and, hence, the renormalization procedure is needed for the boundary-free part only. The 
latter is well-investigated in literature. 

Using formula (3.2), it can be seen that the contribution of the second term on the right of 
Eq. (3.10) to the boundary part of the VEVs vanishes. Prom Eq. (3.8) now one obtains 

where the expressions for the functions F^'^^[g(v)] directly follow from Eqs. (3.9)-(3.11) after the 
integration using formula (3.2): 

F%M\ = (i-2c) [«V'(«') + (i' + j^)«'»(«')9'(<')+ 



i = 0, 1,...,L>-1, (3.13) 
F^'^^igiv)] = F(^^[g{v)]. (3.14) 



By a similar way for the VEVs induced by a single brane in the region z < Za, hy making use 
of expression (2.44) (with replacement zi, — Za), one obtains 

Note that VEVs (3.12), (3.15) depend only on the ratio z/za which is related to the proper 
distance from the plate by equation 

z/za = e'"'^y-''\ (3.16) 



10 



As we see, for the part of the energy-momentum tensor corresponding to the coordinates in the 
hyperplane parallel to the plates one has (r^,)W ~ rini,. Of course, we could expect this result 
from the problem symmetry. It can be seen that the VEVs obtained above obey the continuity 
equation T^j^ = which for the AdS metric takes the form 

^^+1^ {z-^'tE) + DT^ = 0. (3.17) 
In the AdS part without boundary the integration over k can be done using the formula 

k--^,k .°-rf^Vfl--V (3.18) 



Vu^ + A;2 20F V 2 y V 2 

As for the boundary terms, the contribution corresponding to the second term on the right of 
Eq. (3.10) vanishes and one obtains 

iOs\Tt\Os) = ^D^r i^l - -J dnn^-V» [Muz)] , (3.19) 

where the expressions for the functions Z^*-* [9{v)] are obtained from those for F^*) [fl'('i^)]) Eqs. 
(3.13), (3.14), by the replacement v — iv. Now from Eq. (3.13) we have the following relation 
between the vacuum energy density and pressures in the parallel directions (no summation 
over i) 

(0|ro» = (0|T^'|0), i = l,...,D-l. (3.20) 
To evaluate the u-integral in (3.19) we use formula (2.40) and the relations 

+ 2) = - X^^^iD + 2) = - 1 - ? V. - ?1 (3.21) 



A J D + 1 ' ' ' V 2 J \ 2 D + l 



f 

Jo 



( D\ 

dxx^ Mx)Ju-i{x) =\^- -^jMD)- (3.22) 



This yields 



,^ dxx^-'f^^[Mx)] = ^^^, i = 0,l,...,D. (3.23) 



As a result for the boundary-free AdS VEVs one receives 

/oiT^in^ k^-'m' rC-^)r{D/2 + .) 

i^s\Ti \0s) = ^D + mi + u-D/2)- ^^-^^^ 

Hence, the AdS VEVs are proportional to the corresponding metric tensor. Again, we could 
expect this result due to the maximally symmetry of the AdS background. Formula (3.24) can 
be obtained also directly from (2.39), by using the standard relation between the field square 

and the trace of the energy-momentum tensor. 

For a conformally coupled masslcss scalar = 1/2, and by making use of the expressions for 
the modified Bessel functions, it can be seen that {T!')^"-^ = in the region z > Za and 



(4.)-/^r(Z^/2) 7o §^§=^e^ + l 



in the region z < Za- Note that the corresponding energy-momentum tensor for a single Robin 

plate in the Minkowski bulk vanishes [69] and the result for the region z > is obtained by a 
simple conformal transformation from that for the Minkowski case. In the region z < Za this 
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procedure does not work as in the AdS problem one has < z < Za instead of —oo < z < Za 
in the Minkowski problem and, hence, the part of AdS under consideration is not a conformal 
image of the corresponding manifold in the Minkowski spacetime (for a general discussion of 
this question in conformally related problems see Ref. [70]). 

The boundary induced VEVs given by equations (3.12) and (3.15), in general, can not be 
further simplified and need numerical calculations. Relatively simple analytic formulae for the 
branc-induced parts can be obtained in limiting cases. First of all, as a partial check, in the limit 
ko ^ the corresponding formulae for a single plate on the Minkowski bulk are obtained (see 
Ref. [69] for the massless case and Ref. [71] for the massive case). This can be seen noting that in 
this limit v ^ m/ko is large and by introducing the new integration variable in accordance with 
u = vy, we can replace the Bessel modified functions by their uniform asymptotic expansions 
for large values of the order. The Minkowski result is obtained in the leading order. 

In the limit z ^ Za for a fixed kn expressions (3.12) and (3.15) diverge. In accordance with 
(3.16), this corresponds to small proper distances from the brane. The surface divergences in the 
VEVs of the energy-momentum tensor are well-known in quantum field theory with boundaries 
and are investigated for various types of boundary geometries and boundary conditions (see, for 
instance, Refs. [72, 73]). Near the brane the main contribution into the integral over u in Eqs. 
(3.12), (3.15) comes from the large values of u and we can replace the Bessel modified functions 
by their asymptotic expressions for large values of the argument when the order is fixed (see, 
for instance, [67]). To the leading order this yields 

~ ^\-^) 2^7rP+i)/2|l-V^|^+i' (^-'^^ 

(r^)(«) ^ (rO)(«) (i - ^) , (3.27) 

where we use the notation 

Note that the leading terms for the components with z = 0, 1,...,D — 1 arc symmetric with 
respect to the plate, and the g - component has opposite signs for the different sides of the 
plate. Near the brane the vacuum energy densities have opposite signs for the cases of Dirichlet 
{Ba = 0) and non-Dirichlet {Ba 7^ 0) boundary conditions. Recall that for a conformally coupled 
massless scalar the vacuum energy-momentum tensor vanishes in the region z > Za and is given 
by expression (3.25) in the region z < Za- The latter is finite everywhere including the points 
on the plate. 

For the points with the proper distances from the plate much larger compared with the 

AdS curvature radius one has z ^ Za- This limit is important from the point of view of the 
application to the Randall-Sundrum braneworld. Introducing in Eq. (3.12) a new integration 
variable y = zu, by making use of formulae for the Bessel modified functions for small values of 
the argument, and assuming Aa — Bai^ ^ 0, to the leading order we receive 

^^^^ nD/2)T{v)T{v + l)A,-BaA2z) i "^^^ ^ ^^"^^^^ ' ^^"^^^ 

The integral in this formula can be evaluated on the base of formulae 



dxx'^-^K^.ix) = ^r(^ + u)r(--u)r^(-), (3.30) 



/ 

Jo 

00 

,a-l i^'2 



dxx'^-^K^\x) 







r(a) V2 / V2 / V2 
a /a' 
a + 1 yjjo 



o a ( G? 9 
\r 



00 

dxx'^-^Kl{x). (3.31) 
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This yields 



(47r)(^-i)/2 yl„ - S„z/ V2z7 V D + 2I/ + 1 
^ {2v-l)T{D/2 + u + l)T{D/2 + 2v) 

V{v+l)V{D/2 + l/2 + v) ' ^"'•"'^^ 

and the boundary- induced VEVs are exponentially suppressed by the factor exp[2z^fc/5(a — y)]. 
In the limit under consideration the ratio of the energy density and the perpendicular pressure 
is a negative constant. We see that in the case C < Co at large distances from the plate the 
energy density is positive for > \/v and is negative for < 1/u. Combining 

this with the asymptotic behavior near the plate, Eq. (3.26), we conclude that for Ba 7^ and 
\Ba/-^a\ < 1/^^ the energy density is positive near the plate and is negative at large distances 
approaching to zero and, hence, for some value of z it has a minimum with the negative energy 
density. 

Now we turn to the limit z <^ Za for a fixed fc/j. This corresponds to the points near the AdS 
boundary presented by the hyperplane z = with the proper distances from the brane much 
larger compared with the AdS curvature radius. Introducing in Eq. (3.15) a new integration 
variable y = uza and by making use of the formulae for the Bessel modified functions for small 
values of the argument, to the leading order one receives 

/tO\(») fcg^' f^Y^"^ D + 2u-4aD + 2i^ + l) 

^0/ ~ TrD/^T{D/2)\2zaJ r{u)r{u + l) 

X /"dxx^+2-i:^!!M, (3.34) 
Jo i<«)(x) 

(T^)(a) ^(^0°)^"^ z<^za. (3.35) 

In this limit the exponential suppression takes place by the factor exp[(Z) + 2v)kD{y — a)]- The 
ratio of the energy density and the perpendicular pressure is a positive constant. Due to the 
well-known properties of the Bessel modified functions, the integral in Eq. (3.34) is positive for 
small values of the ratio i?„/A„ and is negative for large values and, hence, as in the previous 
case there is a possibility for the change of the sign for the energy density as a function on 
the distance from the boundary. Now by taking into account relation (3.16), from Eqs. (3.32), 
(3.34) we conclude that for large proper distances from the plate to compared with the AdS 
curvature radius the boundary induced parts are exponentially suppressed. 

In the large mass limit, m 3> kn, from (2.11) one has u « m/kD 3> 1. To find the leading 
terms of the corresponding asymptotic expansions for the vacuum energy-momentum tensor 
components we replace the integration variable in Eqs. (3.12) and (3.15), x = vt, and use the 
uniform asymptotic expansions for the Bessel modified functions for large values of the order. 
The main contribution to the integrals over t comes from small values of t. For a fixed z to the 
leading order one receives 

(r^)'"' ~ ^<^>'°'. (3-37) 

and we have an exponential suppression of the vacuum expectation values for large values of the 
mass. 
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For large values of the parameter kf) and fixed proper distances from the plate, k£)\y — o,\ 3> 1, 
from relation (3.16) one has <C z or S> z. Hence, in this limit, which corresponds to strong 
gravitational fields, we have an asymptotic behavior described by Eqs. (3.32), (3.33) in the 
region z > Za and by Eqs. (3.34), (3.35) in the region z < Za with the exponential suppression 
of the boundary induced parts by the factor exp[2i/k£){a — y)] in the first case and by the factor 
exp[(2i^ + D)kD{y — a)] in the second case. 

In figures 1 and 2 we have plotted the vacuum energy density [curves (a)] and ^ - stress 
[curves (b)] induced by a single plate, as functions on koiy — o) for a minimally coupled massless 
scalar field in L> = 3. In figure 1 the cases of Dirichlet (left panel) and Neumann (right panel) 
boundary conditions are presented. The energy density is negative everywhere for Dirichlet 
scalar and is positive for Neumann scalar. It can be seen that for both Dirichlet and Neumann 
cases the boundary-induced energy-momentum tensors violate the strong energy condition. The 
weak energy condition is violated by Dirichlet scalar and is satisfied by Neumann scalar. The 
energy conditions play a key role in the singularity theorems of general relativity and these 
properties may have important consequences in the consideration of the gravitational back re- 
action of vacuum quantum effects. Figure 2 corresponds to the Robin boundary condition with 
Ba/-^a = 1/4 and illustrates the possibility for the change of the sign for the energy density 
as a function on distance from the plate: the energy density is positive near the plate and is 
negative for large distances from the plate. Note that the ratio (T^)(")/(Tq is a coordinate 
dependent function. In accordance with the asymptotic estimates given above, this ratio tends 
to zero for the points near the brane and to the constant values D/(D + 2v) and —D/{2i') in 
the limits z ^ and z <^ Zn. resoectivelv. 
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Figure 1: Boundary induced vacuum densities {Tq)^"''^ jk^'^^ and {T^)^"-^ /k'^'^^ as functions on 
koiy — a) for a minimally coupled massless scalar in D = 3. The curves (a) correspond to the 
energy density and the curves (b) correspond to the ^ - component of the vacuum stress. The 
left panel is for Dirichlet scalar and the right one is for Neumann scalar. 



4 Two-plates geometry 
4.1 Vacuum densities 

In this section we will investigate the VEVs for the energy-momentum tensor in the region 
between two parallel plates. Substituting the corresponding Wightman function from Eq. (2.31) 
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Figure 2: The same as in figure 1 for Robin scalar witli Ba/Aa = 1/4. 
into the mode sum formula (3.5), we obtain 



uD+l Dsk poo 

m^\0) = {Os\Tt\Os) + {Th^^^ - ^D-2/D+i)/2r\l2^) I ("^-^ 



fc V — L 



[UZb, UZ) 



where for a given function g{v) the functions F^'^^[g{v)] are defined in accordance with Eqs. 
(3.9)-(3.11). By using formula (3.2), wc sec that the contribution of the second term on the 
right of Eq. (3.10) vanishes and from Eq. (4.1) one obtains 

UD+I Ds:k poo 

{Q\T^\Q) = {Qs\Tt\Qs) + {Th^""^ - n / duu''-^nau{uZa,uzb)F^'^ G^:^\uZa,uz) , 

2^-i7r2-r(f ) io L J 

(4.2) 

with the functions F'^'\g{v)] from Eqs. (3. 13), (3.14). 

By making use of the Wightman function in form (2.43), we obtain an alternative represen- 
tation of the VEVs: 

(#) 

(4.3) 

This form is obtained from (4.2) by the replacements a ^ h, K^. The vacuum energy- 

momentum tensor in the region between the branes is not symmetric under the interchange of 
indices of the brancs. As it has been mentioned above, the reason for this is that, though the 
background spacetime is homogeneous, due to the non-zero extrinsic curvature tensors for the 
branes, the regions on the left and on the right of the brane are not equivalent. By the same 
way as for the case of a single plate, it can be seen that in the limit ko ^ ^ from formulae (4.2) 
and (4.3) the corresponding results for two plates geometry in the Minkowski bulk (see Refs. 
[69, 71]) arc obtained. 

On the base of formula (4.2), the VEVs in the region Za < z < zi, can be written in the form 
(0|i;*^|0) = (05|7;l05) + {T^t^ + {r^f^ + A{Tt), (4.4) 



h.D+1 D>:k POO 

(0|I^^|0) = {Qs\Th^s) + {T^)^'^ - „^ u ' / duu''-'nUuZa,uzb)F(^^ G^^H 
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with the 'interference' term 



\ ^ I 2D-HD/^T{D/2) Jo 



X \ na,{uZa,uzb)F('^[G<;^HuZa,uz)] - ^^J^ F^"^ [I,{uz)] \ . (4.5) 

[ iy{uzb) J 

Another form wc can obtain by using formulae (3.12) and (4.3). Note that the surface divergences 
are contained in the second and third terms on the right of Eq. (4.4) and expression (4.5) is finite 
for all values Za < z < Z}, (for large values u the subintegrand behaves as i(-Dg2u(za-zi,)) Yov the 
proper distances between the plates much less than the AdS curvature radius, koih — a) <C 1, 
one has (1 — Zajz}^ <^ 1 and the main contribution into the integral in Eq. (4.5) comes from 
large values of u. Replacing the Bessel modified functions by their asymptotic expansions for 
large values of the argument, for the VEVs from Eq. (4.5) to the leading order one receives 

.,rj.D. (i2±l) CnjP + 1) 1 - 2-^-1 [1 - K{Ba)K{B,)] 

^ ^ (47r)(^+i)/2(6-a)«+i k(B„)k(5;,) ^ ^ 



D ' ^' 22^-i7r-D/2r(D/2)(6-a)^+i 



K{Ba)K{Bb)e^ - 1 [ \ Zb- ZaJ \ Zb 



K{Ba) exjp [ t— — —]+ K{Bb)ex.p (t-^ 



(4.7) 



where C,r{z) is the Riemann zeta function and the function k{x) is defined by Eq. (3.28). The 
leading terms given by formulae (4.6) and (4.7) are the same as for the corresponding quantities 
on Minkowski spacetime background. In particular, the 'interference' part of the ^ - component 
does not depend on the curvature coupling parameter. In the limit of large distances between 
the plates, when knih — a) S> 1, introducing a new integration variable v = uzb and expanding 
over Za/zbi to the leading order from Eq. (4.5) one obtains 



2u 



( Z 



A(7;')~4'(5j /'(^j, (4.8) 

where the form of the function g^'^\v) can be found from Eq. (4.5). For v <C 1 one has g(v) ^ 
and from (4.8) it follows that the 'interference' part of the vacuum energy-momentum tensor is 
exponentially suppressed for large interbrane distances. 



4.2 Interaction forces 

Now we turn to the interaction forces between the plates. The vacuum force acting per unit 
surface of the plate at ^; = Zj is determined by the ^ - component of the vacuum energy- 
momentum tensor at this point. The corresponding effective pressures can be presented as a 
sum of two terms: 

P^'^=pf+P%y j = a,b. (4.9) 

The first term on the right is the pressure for a single plate at z = Zj when the second plate is 
absent. This term is divergent due to the surface divergences in the vacuum expectation values. 
The second term on the right of Eq. (4.9), 

P%) = -{tS)^'^-^{tS), z = zj, j,h = a,b, h^j, (4.10) 
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is the pressure induced by the presence of the second plate, and can be termed as an interaction 
force. It is determined by the last terms on the right of formulae (4.2) and (4.3) for the plate at 
z = Za and z = z^ respectively. Substituting into these terms z = zj and using relations 



u 



-B 



J' 



Gl^y{u,u) = A,/u, 



(4.11) 



one has 



P. 



(j) 

(int) 



J,D+1 



2DTrD/2r (f ) 



1-1 rc 



{ [x^ -u'^ + 2m^/kl) B] - D(4C - l)AjBj - A]} 



(4.12) 



Note that due to the asymmetry in the VEV of the energy-momentum tensor, the interaction 
forces acting on the branes are not symmetric under the interchange of the brane indices. By 
taking into account that K^{u)I^{v) — Ky{v)Ii,{u) > and Kl{u)Il{v) — Kl{v)Il{u) < for 
u < V, it can be easily seen that the vacuum effective pressures are negative for Dirichlet scalar 
and for a scalar with Aa = Aj, = and, hence, the corresponding interaction forces are attractive 
for all values of the interplate distance. In figure 3 these forces are plotted as functions on Za/zb 
for massless minimally coupled Dirichlet and Neumann scalars. As it has been shown above, in 
the case of a conformally coupled massless scalar field for the first term on the right of Eq. (4.9) 
one has p^"^ = and p^*^ is determined from Eq. (3.25) with the replacement a ^ b (recall that 
wc consider the region -z^ < z < zj). It can be seen that the corresponding formulae for p^^ 
obtained from Eqs. (4.9), (4.12) coincide with those given in Rcf. [47]. Note that in this case 
p{a) I ^D+i _ pi^) I z^~^^ . Using the Wronskian for the Bessel modified functions, it can be seen 
that 





[Bj{x'^Zj + i/^) - A^] nj^{xza,xzb) = UjZj— In 



li''\xZb)Kj.''\xZa) 



j = a,b, (4.13) 



where = 1, = — 1. This allows us to write the expressions (4.12) for the interaction forces 
per unit surface in another equivalent form: 



(i) 



2^7r^/2r(f)yo 



dx x 



D-l 



l + DB 



2CBj + {l-40Aj 

■^R2(x2z2 + ^2)_^2 



^ 1 

X — In 

OZi 



ii^Hx=a)K::"(xzb) 



'(h) 



1 - 



li''\xZb)KrixZa) 



(»)/ 



(4.14) 



Note that these forces in general are different for j = a and j = b. For Dirichlet scalar the 
second term in the square brackets is zero. 

Let us consider the limiting cases for the interaction forces described by Eq. (4.12). For 
small distances to compared with the AdS curvature radius, koib — a) <^ 1, the leading terms 
are the same as for the plates in the Minkowski bulk: 



4^ = D{1 - 2-^) 



(47r)(^+i)/2(6-a)^+i' 



(4.15) 



in the case of Dirichlet boundary condition on one plate and non-Dirichlet boundary condition 
on the another, and 

P(int) ^(4^)(D+l)/2(^_a)D+l' 



(4.16) 
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Figure 3: The vacuum effective pressures {zh/ Zj)^p'^^^^'^/k'^'^^ , j = a,b, as functions on Za/zi, for 
a massless minimally coupled scalar field in D = 3 determining the interaction forces between 
Dirichlet (left panel) and Neumann (right panel) plates. Curves a and b correspond to the 



pressures on the plates z 



a, b respectively. 



for all other cases. Note that in the first case the interaction forces are repulsive and are 

attractive for the second ease. For large distances between the plates, k^ib — a) ^ 1 (this limit 
is realized in the Randall-Sundrum model), by using the expressions for the modified Bessel 
functions for small values of the argument, one finds 



(a) 
^'(int) 



za\ ^^''^ 4feg+i [(2mV4 - u^)Bl - DjAC - 1)AA - Al] 
2zb J 7r^/2r(Z)/2)r2 (u) [BaV - Aaf 



f 

Jo 



P(int) 



2v 



(4.17) 



Aa + Baiy 



2zJ Aa- Baiy2D-^TTD/2T{D/2)uT^{u) 



dx 



X 



2V+D-1 



[{x" - + 2rr?/kl)Bl - D{AQ - 1)^5;. - J^] . (4.18) 



In dependence of the values for the coefficients in the boundary conditions, these effective pres- 
sures can be either positive or negative, leading to repulsive or attractive forces at large distances. 
To illustrate these possibilities, in figure 4 we have plotted the effective pressures as functions 
on Za/ Zfj for Robin boundary conditions. In the case of the left panel, the vacuum interaction 
forces are attractive for small distances and repulsive for large distances. For the right panel we 
have an opposite situation, the interaction forces are attractive for large distances and repulsive 
for small distances. For the latter case the vacuum interaction forces provide a possibility for 
a stabilization of the interplate distance. The dependence of the vacuum effective pressures on 
the ratio Za/ z^ for various values of the Robin coefficient B^ is given in figures 5 and 6. 



5 Application to the Randall-Sundrum braneworld 

In this section we will consider the application of the results obtained in the previous sections 
to the Randall-Sundrum braneworld model [20] based on the AdS geometry with one extra 
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Figure 4: The same as in figure 3 for mixed boundary conditions. Left panel corresponds to 
Aa = 0, Ba = 1, Ah = 1, Bb = 0.2 and right panel corresponds to Aa = 1, Ba = 0, Ah = 1, 
Bh = 0.2. 
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0.8 0.8 

Figure 5: The D = 3 vacuum effective pressures {zh / Zj)^ p[^^ / k^'^^ , j = a,b, for a massless 
minimally coupled scalar field as functions on Za/zb and Bb- The values for the other Robin 
coefficients are Aa = 1, Ba = 0, and Ah = 1. Left panel corresponds to j = a and right panel 
corresponds to j = b. 
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Figure 6: The same as in figure 5 for the values of Robin coefficients Aa = 0, Ba = 1, and 
Ab = l. 

dimension. The fifth dimension y is compactified on an orbifold, jZ^ of length L, with 
— ^ y ^ L. The orbifold fixed points at y = and y = L are the locations of two 3-branes. 
Below we will allow these submanifolds to have an arbitrary dimension D. The metric in the 
Randall-Sundrum model has the form (2.1) with 

a{y) = kD\y\. (5.1) 

For the corresponding Ricci scalar one has 

R = 4DkD[5{y) - 5{y - L)] - D{D + 1)^. (5.2) 

Let us consider a bulk scalar ip with the curvature coupling parameter C, and with the action 
functional 

S = \j d'^xdyy/lil [g^'^di^dkip - [m^ + ci5(y) + C2(5(y - L) + CR] if^] , (5.3) 

where m is the bulk mass, ci and C2 are the brane mass terms on the branes y = and y = L 
respectively. The field equation obtained from (5.3) is 

V,VV + K + ^^^(y) + C2<5(y -L) + CR]ip = 0. (5.4) 

The corresponding eigenfunctions can be written in form (2.7) and (2.8), where fniv) is a solution 
to the equation 



D(T_^ f -Dadfn 

ay 



^) + [rr? + {cl+ADCkD)5{y) + {c2-^DCkD)5{y-L)- 
-D{D + l)Ckl] U{y) = my^Uy). (5.5) 



To obtain the boundary conditions for the function fn{y) we integrate (5.5) first about y = 
and then about y = L. Assuming that the function /n(y) is continuous at these points, we 
receive 



hm^ 

dy 



lZ%^['}[y, = [cj - {-lymko] fniyj), j = 1, 2, yi = 0, y2 = L , (5.6) 
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with e > 0. First consider the case of the untwisted scalar field for which fn{y) = fn{—y)^ 
and the solution to Eq. (5.5) for y 7^ 0, L is given by expression (2.10) with z = e'^^I^L The 
boundary conditions which follow from relations (5.6) are in form (2.4) with (see also Ref. [43] 
for the case ci = C2 = 0) 



^ = -bc^ + 4DCkD), ^ = -l{-C2 + mkD), (5.7) 
i>o ^ rib ^ 



and respectively 



^ = ^ [D(l - 40 - ci/M , ^ = ^[^(1-4C)+C2/M- (5.8) 

For the twisted scalar fn{—y) = —fn{y) and from the Z2 identification on one has /n(0) = 
fn{L) = and, hence, Ba = Bf, = 0, which correspond to Dirichlet boundary conditions on both 
branes. The normalization condition for the eigenfunctions /n(y) now has the form 



dye^^-^>^y^fn{y)fn'iy) = Snn'- (5.9) 

L 



As a result the normalization coefficient c„ differs from (2.15) by additional factor 1/2. Hence, 
the Wightman function in the Randall-Sundrum braneworld is given by formulae (2.31) or 
equivalently by (2.43) with an additional factor 1/2 and with Robin coefficients given by Eq. 
(5.7). Similarly, the VEV of the energy-momentum tensor is defined by formulae (4.2) or (4.3) 
by an additional factor 1/2. Note that global quantities, such as the total Casimir energy, are the 
same. Motivated by the possibility for the stabilization of extra dimensions by quantum effects, 
the one-loop Casimir energy of various braneworld compactifications has been investigated by 
several authors for both scalar and fermion fields (see references cited in Introduction). It has 
been shown that for the simultaneous solution of the stabilization and hierarchy problems a fine 
tuning of the model parameters is essential. From the point of view of backrcaction of quantum 
effects the investigation of local densities is also of considerable interest. For a conformally 
coupled massless scalar field the corresponding results can be obtained from the results in the 
Minkowski spacetime by using the standard formula for conformally related problems. For the 
case of plane-symmetric line element (2.1) with an arbitrary function (T(y) this has been done in 
Ref. [47] by using the results from [69]. Recently the energy-momentum tensor in the Randall- 
Sundrum braneworld for a bulk scalar with zero brane mass terms ci and C2 is considered in Ref. 
[49] . This paper appeared when our calculations were in progress. Note that in [49] only a general 
formula is given for the unrenormalized VEV in terms of the differential operator acting on the 
Green function. In our approach the application of the generalized Abel-Plana formula allowed 
us to extract manifestly the part due to the AdS bulk without boundaries and for the points 
away from the boundaries the renormalization procedure is the same as for the boundary-free 
parts. The latter is well-investigated in literature. In addition, the boundary induced parts are 
presented in terms of exponentially convergent integrals convenient for numerical calculations. 



6 Conclusion 

The natural appearance of AdS in a variety of situations has stimulated considerable interest 
in the behavior of quantum fields propagating in this background. In the present paper we 
have investigated the Wightman function and the vacuum expectation value of the energy- 
momentum tensor for a scalar field with an arbitrary curvature coupling parameter satisfying 
Robin boundary conditions on two parallel plates in AdS spacetime. The application of the 
generalized Abel-Plana formula to the mode sum over zeros of the corresponding combinations 
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of the cylinder functions allowed us to extract the boundary-free AdS part and to present the 
boundary induced parts in terms of exponentially convergent integrals. In the region between 
two plates the Wightman function is presented in two equivalent forms given by Eqs. (2.31) 
and (2.43). The first terms on the right of these formulae are the Wightman function for AdS 
bulk without boundaries. The second ones are induced by a single plate and the third terms 
are due to the presence of the second plate. The expectation values for the energy-momentum 
tensor are obtained by applying on the Wightman function a certain second order differential 
operator and taking the coincidence limit. For the case of a single plate geometry this leads to 
formula (3.12) for the region z > Za and to formula (3.15) for the region z < Za- As we could 
expect from the problem symmetry the part of this tensor corresponding to the components 
on the hypcrplanc parallel to the plate is proportional to the corresponding metric tensor. On 
the boundary the vacuum energy-momentum tensor diverges, except the case of a conformally 
coupled massless scalar. The leading terms of the corresponding asymptotic expansion near the 
boundary are given by expressions (3.26) and (3.27). These terms are the same as for a plate 
in the Minkowski bulk. They do not depend on the Robin coefficient, have different signs for 
Dirichlet and Neumann scalars and vanish for a conformally coupled scalar. The coefficients for 
the subleading asymptotic terms will depend on the AdS curvature radius, Robin coefficient and 
on the mass of the field. For large proper distances from the plate to compared with the AdS 
curvature radius, kolu — ci\ S> 1, the boundary induced energy-momentum tensor vanishes as 
exp[2iyku{a — y)] in the region y > a and as cxp[k£){2i' + D)(y — a)] in the region y < a. The same 
behavior takes place for a fixed y—a and large values of the parameter ko- In the large mass limit, 
m S> ko, the boundary parts are exponentially suppressed. Note that here we consider the bulk 
energy-momentum tensor. For a scalar field on manifolds with boundaries in addition to this 
part, the energy-momentum tensor contains a contribution located on on the boundary (for the 
expression of the surface energy-momentum tensor in the case of arbitrary bulk and boundary 
geometries see Ref. [74]). As it has been discussed in Refs. [73, 69, 65, 75, 76, 74], the surface 
part of the energy-momentum tensor is essential in considerations of the relation between local 
and global characteristics in the Casimir effect. The vacuum expectation value of the surface 
energy-momentum tensor for the geometry of two parallel branes in AdS bulk is evaluated in 
Ref. [77] . It is shown that for large distances between the branes the induced surface densities 
give rise to an exponentially suppressed cosmological constant on the brane. In particular, in 
the Randall-Sundrum model the cosmological constant generated on the visible brane is of the 
right order of magnitude with the value suggested by the cosmological observations. 

For the case of two-plates geometry the vacuum expectation value of the bulk energy- 
momentum tensor is presented as a sum of purely AdS, single plates, and 'interference' parts. 
The latter is given by Eq. (4.5) and is finite for all values Za < z < zi,. In the limit Za Zh 
the standard result for two parallel plates in the Minkowski bulk is obtained. For two-plates 
case the vacuum forces acting on boundaries contain two terms. The first ones are the forces 
acting on a single boundary when the second boundary is absent. Due to the well-known surface 
divergencies in the expectation values of the energy-momentum tensor these forces are infinite 
and need an additional regularization. The another terms in the vacuum forces are finite and 
are induced by the presence of the second boundary and correspond to the interaction forces 
between the plates. These forces are given by formula (4.12) with j = a,b for the plate at z = Za 
and z = Zf, respectively. For Dirichlet scalar they are always attractive. In the case of mixed 
boundary conditions the interaction forces can be either repulsive or attractive. Moreover, there 
is a region in the space of Robin parameters in which the interaction forces are repulsive for 
small distances and are attractive for large distances. This provides a possibility to stabilize in- 
terplate distance by using the vacuum forces. For large distances between the plates, the vacuum 
interaction forces per unit surface are exponentially suppressed by the factor exp[2z/fcx)(a — b)] 
for the plate at y = a and by the factor exp[A;£)(2i^ + D){a — b)] for the plate at y = b. In Section 
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5 we give an application of our results for two plates case to the Randall-Sundrum braneworld 
with arbitrary mass terms on the branes. For the untwisted scalar the Robin coefficients are 
expressed through these mass terms and the curvature coupling parameter. For the twisted 
scalar Dirichlet boundary conditions are obtained on both brancs. Note that in this paper we 
have considered boundary induced vacuum densities which are finite away from the boundaries. 
As it has been mentioned in Ref. [78], the same results will be obtained in the model where 
instead of externally imposed boundary condition the fluctuating field is coupled to a smooth 
background potential that implements the boundary condition in a certain limit [79]. 
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